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I.  SUMMARY 

This  report  describes  the  operation,  maintenance  and  research  activities  at 
the  Norwegian  Seismic  Array  (NORSAR)  for  the  period  1  April  -  30  September 
1984. 

The  uptime  of  the  NORSAR  online  detection  processor  system  has  averaged 
95.1  per  cent  as  compared  to  95.6  for  the  previous  period.  Almost  half  the 
downtime  came  in  July  and  was  caused  by  power  supply  failure  on  the  MODCOMP 
computer.  A  total  of  2160  events  were  reported  in  this  period,  giving  a 
daily  average  of  11.8  events.  Frequent  tests  and  interruptions  of  com¬ 
munications  systems  for  02C  and  03C  resulted  in  high  error  figures  during 
week  15.  These  tests  were  connected  with  the  malfunction  in  the  MODCOMP 
computer,  caused  by  a  move  in  connection  with  a  reorganization  of  the  com¬ 
puter  hall. 

As  a  result  of  unsatisfactory  performance  of  the  old  alarm  system,  a  new 
alarm  has  been  installed  in  the  computer  room.  The  new  system  has  resulted 
in  safer  reporting  of  breakdowns  in  the  online  system  and  will  hopefully 
result  in  an  improved  uptime.  A  new  telex  connection  has  been  established 
enabling  an  automatic  transmission  of  bulletin  data  to  AFTAC  3  times  a  day, 
7  days  a  week. 

Section  V  describes  maintenance  activities  in  the  field  and  at  the  NORSAR 
Maintenance  Center. 

The  research  activity  is  briefly  described  in  section  VII.  Section  VII. 1 
discusses  source  solutions  from  long-period  waveform  inversion.  Section 
VII. 2  presents  a  study  on  path  effects  and  long-period  station  residuals  in 
the  GDSN  network.  Section  VII. 3  presents  a  preliminary  report  on  a  method 
of  dynamic  ray  tracing,  while  VII. 4  describes  regional  arrays  and  optimum 
data  processing  schemes.  Section  VII. 5  takes  up  high-frequency  P  wave  atte¬ 
nuation  from  Central  Asia  to  Norway.  In  Section  VII. 6  some  initial  results 
are  summarized  with  regard  to  possible  improvements  in  NORSAR  on-line  event 
detection  performance.  Section  VII. 7  contains  a  review  of  recent  develop¬ 
ments  concerning  the  new  NORESS  small-aperture  array. 


II.  OPERATION  OF  ALL  SYSTEMS 

II. 1 _ Detection  Processor  (DP)  Ope ration 

There  have  been  126  breaks  in  the  otherwise  continuous  operation 
of  the  NORSAR  online  system  within  the  current  6-month  reporting 
interval.  Almost  half  the  downtime  came  in  July  and  was  caused  by 
power  trouble  on  the  MODCOMP  computer.  The  uptime  percentage  for 
the  period  is  95.1  as  compared  to  95.6  for  the  previous  period. 

Fig.  II. 1.1  and  the  accompanying  Table  II. 1.1  both  show  the  daily 
DP  downtime  for  the  days  between  1  April  and  30  September  1984. 
The  monthly  recording  times  and  percentages  are  given  in  Table 
II. 1.2. 

The  breaks  can  be  grouped  as  follows: 


a)  MODCOMP  failure  47 

b)  Stops  related  to  program  work  or  error  8 

c)  Maintenance  stops  8 

d)  Power  jumps  and  breaks  5 

e)  Stops  related  to  system  operation  2 

f)  TOD  error  correction  39 

g)  Communication  lines  17 


The  total  downtime  for  the  period  was  215  hours  and  1  minute.  The 
mean-time-between-failures  (MTBF)  was  1.4  days,  the  same  as  for 
the  previous  period. 
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Month 

DP  Uptime 
(hrs) 

DP  Uptime 
(X) 

No.  of 

DP  breaks 

No.  of  days 
with  breaks 

DP 

MTBF* 

(days) 

Apr 

714.08 

99.2 

16 

9 

1.8 

May 

743.48 

99.9 

10 

9 

2.8 

Jun 

715.88 

99.4 

17 

11 

1.7 

Jul 

637.63 

85.7 

43 

20 

0.6 

Aug 

716.98 

96.4 

10 

10 

2.7 

Sep 

649.05 

90.1 

30 

18 

0.9 

4177.10 

95.1 

126 

77 

1.4 

Mean-time-between-f allures  »  (total  uptioe/no.  of  up  Intervals) 


Table  II. 1.2 


Online  system  performance 


1  April 


30  September  1984 


II. 2 _ Array  communications 

Table  II. 2. I  reflects  the  performance  of  the  communications 
systems  throughout  the  reporting  period.  Numbers  marked  with  asterisks, 
indicating  high  error  rates,  are  commented  on  in  the  following 
paragraphs . 

Modcomp  stops  reported  in  the  previous  summary  continued  in  April. 

Until  two  communications  interface  cards  were  swapped  in  Modcomp  (14 
April),  02C  and  03C  communications  systems  were  suspected  as  sources 
of  the  problems.  Frequent  tests  and  interruptions  of  these  com¬ 
munications  systems  resulted  in  high  error  figures  during  week  15. 

On  July  10  the  MODCOMP  processor  was  moved  to  another  location  in 
connection  with  preparations  for  installing  new  equipment  in  the  com¬ 
puter  hall.  After  reconnection  and  restart  the  Modcomp  stopped  fre¬ 
quently  as  a  result  of  power  supply  failure.  In  addition  a  bad  contact 
in  the  plug  between  the  high  speed  modem  and  the  Modcomp  was  found. 

Although  the  performance  has  since  improved,  we  have  temporarily  lost 
the  ability  to  automaticallyinitiate  resynchronization  of  the  lines 
after  restart.  This  is  not  of  vital  importance,  as  resynchronization 
may  also  be  initiated  "manually"  by  the  operator  using  a  special 
option  in  the  online  program.  The  problem  will  be  taken  care  of  when 
other  high  priority  tasks  have  been  finished  (such  as  modification  of 
the  1950  Interface  Handler  between  the  Modcomp  and  the  IBM  4341, 
etc. ) . 
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Summary  by  month: 

Apr  Problem  with  Modcomp  stops  continued.  02C,  03C  com¬ 

munications  systems  were  suspected.  In  connection  with 
attempts  to  solve  the  problem  the  systems  were  fre¬ 
quently  affected  (tests,  etc.),  causing  discontinuity  in 
their  operation.  Consequently,  high  error  figures 
(week  15),  21.4%. 

May  Two  communication  interface  cards  were  swapped  in  Mod¬ 

comp  (14  Apr)  improving  the  performance  drastically. 
Hence  02C,  03C  back  to  normal  operation.  02B  was  down 
between  week  3  and  4  due  to  lack  of  power  to  the 
subarray. 

June  As  a  result  of  power  loss  (after  lightning)  04C  was  down 
between  5  and  25  June. 

July  In  the  beginning  of  July  data  equipment,  including  the 
Modcomp,  were  moved  to  other  locations  in  the  computer 
hall,  causing  Modcomp  problems  with  many  stops  and  high 
error  figures,  as  indicated  In  Table  II.2.1. 

Aug  04C.  Error  in  data  stream  last  weeks,  average  0.13%. 

NTA  Lillestr0m  and  Hamar  checked  line.  Irregularities 
were  not  discovered.  Still  minor  problems  with  the 
Modcomp . 

Sept  04C  was  visited  (Oct  3)  as  communications  system  still 
erroneous.  Modem  replaced  without  improving  the 
performance.  Test  indicated  line  problems. 
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Also  02B  communications  system  is  erroneous.  Test  result 
gave  us  reason  to  believe  the  path  towards  the  CTV 
caused  data  errors. 

02C  communications  system  was  down  between  11  and  20 
September  after  cable  rerouting  and  movement  of  equip¬ 
ment  caused  problems . 

Between  15  and  17  September  03C  communications  system 
was  switched  out  due  to  resynch,  problems. 

06C  (NORSAR  subarray)  will  be  made  operational  again  as 
soon  as  the  line  and  data  communications  equipoment  are 
checked . 

Miscellaneous 

Operation  of  the  analog  lines  from  NORESS  stations  was  terminated  in 
the  beginning  of  May.  The  only  analog  line  in  operation  is  the  one 
from  the  02B  (Telemetry  station)  whose  performance  has  been  satisfac¬ 
tory. 


0. A.  Hansen 


Sub-  Apr  84  (4)  May  84  (5)  Jun  84  (4)  Jul  84  (4)  Aug  84  (5)  Sep  84  (4)  Average 
array  (2-29.4)  (30.4-3.6)  (4. 6-1. 7)  (2-29.7)  (30.7-2.9)  (3-30.9)  i  yr 
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III.  ARRAY  PERFORMANCE 

III.l  Event  Processor  Operation 

In  Table  III.l  some  monthly  statistics  of  the  Event  Processor 
operation  are  given: 


Teleseismic 

Core  Phases 

Sum 

Dally 

Apr  84 

295 

58 

353 

11.8 

May  84 

321 

78 

399 

12.9 

Jun  84 

275 

52 

327 

10.9 

Jul  84 

308 

46 

354 

11.4 

Aug  84 

331 

51 

382 

12.3 

Sep  84 

256 

89 

345 

11.5 

1786 

374 

2160 

11.8 

Table  III.l 

Event  Processor  statistics.  Apr  -  Sep  1984 


B.Kr.  Hokland 
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IV.  IMPROVEMENTS  AND  MODIFICATIONS 

IV. 1  NORSAR  on-line  system  using  4331/4341  and  MODCOMP  Classic 
We  refer  to  the  detection  processor  operation  statistics  for 
detailed  Information  about  uptime,  communications  and  processing 
performance. 

During  this  reporting  period  we  had  one  serious  breakdown  in  the 
on-line  system  caused  by  the  MODCOMP.  The  breakdown  occurred 
after  we  moved  the  computer  as  a  result  of  a  reorganization  of 
the  computer  hall.  Service  people  from  the  computer  company  used 
two  days  to  apparently  solve  the  problem.  However,  we  soon  saw  a 
lot  of  spikes  In  the  data,  and  a  few  days  later  the  main  problem 
was  located  in  the  power  supply.  After  a  new  power  supply  was 
installed  the  system  has  been  running  very  well.  A  previous 
problem  with  0.5  sec  data  drop  (mentioned  In  the  previous  report) 
also  disappeared  after  the  replacement  of  the  power  supply. 

As  a  result  of  bad  performance  of  the  old  alarm  system  we  have 
installed  a  new  alarm.  The  system  calls  the  operator  on  duty  when 
there  is  a  power  break  or  stop  in  the  communications  system.  The 
new  alarm  system  has  resulted  in  a  safer  reporting  of  breakdowns 
in  the  on-line  system  and  will  hopefully  result  in  an  improved 
upt lme . 

A  new  telex  connection  is  now  available  through  our  SERIES/1 
computer.  This  connection  makes  it  possible  for  us  to  create  a 
telex  on  the  computer  and  then  transmit  the  telex  directly  into 
the  telex  network.  As  a  result  of  an  AFTAC  request  for  more  fre¬ 
quent  transmission  of  bulletins,  we  have  Implemented  an  automatic 
procedure  for  DPX  and  EPX  transmissions,  using  the  new  telex 
Interface.  During  the  past  two  months  our  new  telex  system  with 


automatic  transmission  of  bulletins  has  been  tested  with 
transmission  3  times  a  day  and  7  days  a  week.  The  data  are 
prepared  by  the  event  processing  system  and  up-to-date  bulletins 
are  transmitted  at  0100,  0900  and  1700  GMT. 

The  NORSAR  on-line  system  is  now  completely  automatic  after  this 
telex  installation  with  automatic  data  collection,  real-time  data 
processing,  plotting  of  events  and  transmission  of  bulletins. 

R.  Paulsen 


In  the  period  P.W.  Larsen  has  devoted  much  time  to  tasks  related 
to  NORESS.  He  has  participated  in  meetings  with  consultants  and 
contractors,  field  surveys  and  field  installations.  Also  admi¬ 
nistration  and  support  in  connection  with  the  U.S.  team  from  San- 
dia  and  Comsat,  etc.,  has  taken  part  of  his  time. 


Visits  to  the  subarrays  have  been  limited  to  cases  of  more  or 
less  "acute"  character,  affecting  the  operation  of  a  whole 
subarray,  such  as  power  loss,  failing  modems,  or  other  situations 
and  incidents  which  needed  NMC  assistance.  On  the  other  hand,  per 
sonnel  at  the  computer  center  have  been  able  to  adjust  most  LP 
seismometers  which  were  out  of  adjustment .  By  means  of  data 
collection  and  analysis  programs  reduced  performance  and  malfunc¬ 
tions  could  be  discovered. 


Table  V.l  indicates  the  visits  made  to  subarrays  in  the  period. 


I 


Subarray/ 


Task 


Date 


area 


V.2 _ Improvements  and  modifications 

No  changes  have  been  made  In  the  NORSAR  or  NORESS  arrays. 

Recording  of  analog  data  from  the  provisional  NORESS  array  was 
terminated  at  the  beginning  of  May  1984. 

V.3 _ Array  status 

As  of  30  September  1984  the  following  channels  deviated  from 
tolerance: 

01B  05,  08,  06C  (not  in  operation) 

03C  08 

04C  09 

01A  01  8  Hz  filter 

02 

04  AH.  30  dB 


O.A.  Hansen 
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VI.  DOCUMENTATION  DEVELOPED 

Doornbos,  D.J.:  Solutions  and  station  residuals  from  long-period 

waveform  Inversion  of  deep  events.  Submitted  for  publication. 

Husebye,  E.S.,  S.F.  Ingate  &  E.  Thoresen:  Seismic  arrays  for  every¬ 
one.  Submitted  for  publication. 

Husebye,  E.S.  and  S.F.  Ingate:  Seismic  arrays  -  a  new  rennaissance. 
Submitted  for  publication. 

Ingate,  S.F.,  E.S.  Husebye  and  A.  Christoff ersson:  Regional  arrays 
and  optimum  processing  schemes.  Submitted  for  publication. 

Tronrud,  L.B.:  Semiannual  Tech.  Summary,  1  Oct  1983  -  31  Mar  1984, 
NTNF/NORSAR,  Kjeller,  Norway. 


L.B.  Lough ran 
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VII.  SUMMARY  OF  TECHNICAL  REPORTS/PAPERS  PREPARED 

VII.  1 _ Source  solutions  from  long-period  waveform  Inversion 

Seismic  waveform  inversion  techniques  are  being  applied  to  deter¬ 
mine  source  parameters,  and  to  determine  the  parameters  describ¬ 
ing  earth  structure.  This  requires  (1)  a  suitable  source  para¬ 
meterization,  and  (2)  a  suitable  device  to  account  for  earth 
structure  effects  in  waveform  modelling.  Moment  tensor  techniques 
offer  a  partial  solution  to  the  first  problem  since  the  source 
type  need  not  be  specified  a  priori,  whereas  model-dependent 
constraints  can  still  be  used  in  the  inversion.  A  remaining 
problem  is  that  a  moment  tensor  method  in  practice  Involves  a 
point  source  inversion  technique,  i.e.,  solutions  are  obtained 
with  a  reference  to  a  point  (^o,t0).  In  previous  work  aimed  at 
the  determination  of  source  finiteness,  a  particular  par- 
torization  of  the  moment  tensor  stress  glut  was  used: 

mjkd.O  -  Mjk  f(£,-t)  (1) 

and  the  source  finiteness  is  contained  in  f(^,T>.  However,  the 
effect  of  source  finiteness  can  usually  be  neglected  in  long- 
period  waves  from  sources  of  up  to  moderate  size.  Alternatively, 
if  the  effect  is  known  (or  assumed),  it  can  be  directly  included 
in  the  Green's  functions.  In  these  circumstances,  a  useful  ap¬ 
proximation  to  the  far-field  response  is  a  sum  of  terms  of  the 
form 


i 

UiCx.t)  -  Mjk  Gjjk  (^j.x.t-TQ-STi)  (2) 

where  Gjfk*  is  a  modified  Green's  function,  and  the  time  delay 
depends  on  the  mislocation  of  the  source.  In  practice,  the 
use  of  an  imperfect  earth  model  in  computing  Gj>k*  also  leads  to 
a  time  residual,  so  we  can  write 
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5xi  -  CiT  F(1)  +  di  (3) 

where  £(i)  is  the  mislocation  vector,  ^  is  a  slowness  vector, 
and  di  is  often  called  a  station  residual. 

The  straightforward  approach  to  solve  for  6x£  in  equation  (2)  is 
by  a  correlation  method,  e.g.,  Ward  (1983).  In  the  iterative 
method  we  use,  Mj^  is  estimated  by  a  least  squares  inversion  of 
the  normalized  records;  this  estimate  is  used  in  forming  syn¬ 
thetics  which  are  correlated  with  the  observed  waveforms  to  de¬ 
termine  5ti ,  etc*  With  shallow  sources,  where  the  observed 
records  are  modelled  by  a  superposition  of  direct  and  depth 
phases,  the  correlation  method  should  be  applied  for  a  range  of 
depths. 

Fig.  VII. 1.1  illustrates  the  applicaiton  of  this  method  to  a 
deep  event  in  the  Fiji  Islands  region.  Rather  than  showing  the 
complete  moment  tensor,  the  equal  area  projection  to  the  right 
shows  only  a  double  couple  approximation.  However,  results  on 
deviations  from  the  double  couple  for  12  events  are  summarized 
in  Fig.  VII. 1.2.  Figs.  VII. 1.2a,  b  and  c  were  obtained  by 
slightly  different  Inversion  procedures,  and  the  result  in  some¬ 
what  different  deviations.  Our  presumably  best  solutions  are  in 
Fig.  VII. 1.2c,  and  we  conclude  that  at  this  stage  there  is  no 
compelling  evidence  for  significant  deviations  from  the  double 
couple  mechanism  of  these  deep  earthquakes  (c.f.  Giardinl,  1984). 

The  method  also  leads  to  an  estimate  of  the  centroid  time  shift 
Ax0  which  is  one  of  the  components  of  the  mislocation  vector,  and 
which  is  related  to  source  flnlteness.  The  centroid  time  shifts 
are  shown  in  Fig.  VII. 1.3  together  with  the  time  difference  be¬ 
tween  the  short-period  P  onset  and  the  long-period  P  correlation 
lag.  The  theoretical  relation  would  be  of  the  form 


6x0  -  b(M/Ao)1/3 


(4) 


and  for  a  circular  fault  with  directivity  coefficient  d  ■  1.5 
(Doornbos,  1984)  and  relevant  P  and  S  velocities:  b  ■  17.2x10-9 
where  M  is  in  dyne. cm  and  Ax  in  bar.  Comparison  with  the  results 
in  Fig.  VII. 1.3  then  suggests  that  the  data  are  consistent  with 
a  stress  drop  in  the  range  10-100  bar. 


D.J.  Doornbos 
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Fig.  VII. 1.1  Vertical  and  transverse  component  records  from  sta¬ 
tions  of  the  GDSN  network,  showing  P,  pP,  SH  and 
sSH  from  a  deep  m^  5.6  event  in  the  Fiji  Islands 
region.  Synthetic  record  sections  for  the  moment 
tensor  point  source  solution  are  plotted  over  the 
observed  records.  To  the  right,  the  equal  area  pro¬ 
jection  of  a  double  couple  approximation  to  the 
moment  tensor. 
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Fig.  VII. 1.2  Deviation  from  the  double  couple  versus  the  seismic 
moment  for  12  deep  events  in  the  Fiji  Islands  re¬ 
gion.  The  deviation  is  given  by  X.3  |  |  where 

|  |  and  |  X3  |  are  the  largest  and  smallest  absolute 

eigenvalue  of  the  moment  tensor.  A  positive/negative 
deviation  corresponds  with  being  the  compress ional/ 
tensional  eigenvalue.  The  results  in  a,  b  and  c  repre¬ 
sent  different  moment  tensor  inversions  but  the  best 
solutions,  with  standard  deviations,  are  given  in  c. 


MOMENT  (dyne  cm) 


Fig.  VII. 1.3  0:  Centroid  time  shift  versus  moment  of  the  deep 

events  referred  to  In  Fig.  VII. 1.2.  t:  Averaged 
difference  between  short-period  P  onset  and  long- 
period  P  correlation  lag.  The  solid  line  is  a  theo 
retlcal  relation  of  the  form  Ax0  ■  c  with 

c  -  3.6x10~9  s/(dyne.cra)l/3. 


VII. 2 


Some  path  effects  and  long-period  station  residuals 
In  the  GDSN  network 
In  present-day  routine  methods  of  source  analysis  (Dziewonski  and 
Woodhouse,  1983),  intrinsic  source  parameters  and  the  mislocatlon 
of  the  source  can  be  determined  simultaneously  from  long-period 
records,  for  a  fixed  earth  model.  It  has  been  argued  that  lateral 
relocation  of  events  is  often  needed  not  so  much  for  the  purpose 
of  correcting  real  mislocations ,  but  rather  as  a  device  to  com¬ 
pensate  for  the  travel  time  or  phase  delay  effects  of  anomalous 
earth  structure.  The  implication  of  this  argument  is  in  accor¬ 
dance  with  other  findings  indicating  that,  away  from  critical 
zones,  the  main  effect  of  velocity  anomalies  is  often  on  the  tra¬ 
vel  times  of  ray  contributions  to  the  long-period  body  waves 
(e.g.,  Wallace,  1983).  A  logical  consequence  of  this  line  of 
reasoning  is  to  try  and  simultaneously  determine  intrinsic  source 
parameters  and  travel  time  residuals  rather  than  source  location. 
The  formulation  of  this  procedure  is  outlined  in  a  comparison 
contribution  (see  also  Doornbos,  1984).  We  analyze  the  travel 
time  residuals  by  writing  them  in  the  form 

=  iijT  I(l)ij  +  di  (1> 

where  F.(l)  is  the  mislocatlon  vector  of  the  source,  C  is  a 
slowness  vector,  and  i  and  j  denote  the  source  and  receiver, 
respectively.  Our  present  analysis  of  station  residuals  uses  PREM 
as  the  reference  model.  Since  we  keep  the  velocity  model  fixed, 
all  propagation  path  effects  are  lumped  into  the  station  term  d^. 
The  station  term  is  different  for  P  and  S  waves,  but  for  the 
moment  we  will  assume  that  it  is  otherwise  constant. 

We  can  eliminate  the  depth  mislocatlon  term  in  equation  (1). 
Furthermore,  in  the  present  work  we  have  kept  the  lateral  loca¬ 
tion  of  the  source  fixed.  The  only  remaining  term  is  then  the 
centroid  time  shift  At0)j  -  (t0-T0)j.  There  is  a  natural  nonuni¬ 
queness  of  the  division  of  time  delay  in  terms  in  equation  (1). 
In  terms  of  our  simplified  relation,  it  means  that  the  only 
constraint  on  the  average  of  station  residuals  d  is  that 


Using  a  data  set  from  12  deep  events  in  the  Fiji  Islands  region, 
we  have  measured  the  time  difference  between  the  short-period  P 
onset  and  the  long-period  correlation  lag.  The  average  of  these 
differences  for  each  event  is  taken  to  be  a  measure  of  AT0>j;  the 
average  over  all  events  is  then  a  measure  of  At0.  The  system  of 
equations  (1)  has  been  solved  with  this  constraint  for  Atq,  and 
the  resulting  station  residuals  for  P  and  S  waves  are  given  as  a 
function  of  epicentral  distance  in  Fig.  VII.2.1.  Typical  obser¬ 
vational  features  are  the  base  line  effect  for  P  and  S,  the  trend 
with  epicentral  distance,  especially  for  S,  and  the  travel  time 
anomalies  for  S  at  some  of  the  stations.  All  of  these  features 
need  further  investigation,  but  modifications  to  the  reference 
model  (including  the  absorption  band  used)  seem  inevitable. 


We  note  that  not  only  the  time  residuals,  but  also  the  amplitude 
residuals  are  significant.  The  amplitude  residue  is  defined  as 
the  amplitude  ratio  observed/predicted  where  the  prediction  Is 
based  on  the  moment  tensor  solution  and  the  residence  is  an 
average  over  the  events,  i.e.,  it  is  the  least  squares  estimate 
based  on  the  individual  residuals  as  obtained  from  each  moment 
tensor  solution.  Fig.  VII. 2. 2  shows  the  amplitude  residues  of  P 
and  Fig.  VII. 2. 3  those  of  SH,  but  we  have  plotted  the  residues 
for  direct  (P  and  SH)  and  depth  phases  (pP  and  sSH)  separately. 
There  is  no  clear  correlation  with  distance,  but  there  is  a  clear 
difference  between  the  residuals  for  direct  and  depth  phases, 
especially  for  SH.  In  light  of  these  results,  we  would  have  to 
consider  necessary  modifications  of  in  particular  the  upper 
mantle  in  the  source  region. 


D. J.Doornbos 
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Fig.  VII. 2. 2  Amplitude  residues  for  long-period  P  (•)  and  pP 
(o)  at  the  GDSN  stations  as  in  Fig.  VII. 2.1.  The 
residue  is  defined  as  the  amplitude  ratio  observed/ 
predicted.  Stations  are  positioned  at  their  average 
distance  from  the  events. 
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Fig.  VII. 2. 3  Amplitude  residues  for  long-period  SH  (•)  and 
(o)  at  the  GDSN  stations  as  in  Fig.  VII.2.1.  0 
details  as  in  Fig.  VII. 2. 2. 
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VII. 3  Preliminary  report  on  a  method  of  dynamic  ray  tracing 
The  project,  started  in  October  1983,  is  still  in  progress.  The 
objective  is  numerical  Implementation  of  the  methods  of  dynamic 
ray  tracing  (of  seismic  waves)  developed  by  the  author.  These 
methods  encompass  (1)  propagation  through  caustic  regions;  (2) 
edge  diffraction;  (3)  the  effects  of  anisotropy  and  prestress. 

Algorithms  and  FORTRAN  programs  for  kinematic  ray  tracing  and 
dynamic  ray  tracing  have  been  worked  out  by  A.  Hanyga.  Prelimi¬ 
nary  testing  against  simple  analytical  models  (SATM2,  SATM3)  has 
been  carried  out  by  A.  Hanyga  on  IBM  at  NORSAR,  partly  in  colla¬ 
boration  with  O.A.  Sandvin.  Testing  against  realistic  models  has 
been  done  by  Dr.  Jan  Pajchel  at  the  Seismological  Observatory, 
University  of  Bergen,  on  VAX. 

The  programs  trace  rays  and  signals  in  arbitrary  user-supplied 
models.  The  models  are  represented  by  model  data  files  and  by  a 
library  of  subprograms  called  by  the  ray  tracing  programs.  We  now 
describe  various  aspects  of  the  methods  applied. 

Dynamic  ra£  .tracing  thorough  caustic^  re^ion^ 

It  has  been  shown  by  Hanyga  (1984,  1984a,  1984b,  1984c)  that 
dynamic  ray  tracing  can  be  generalized  to  cover  the  case  of  WKBJ 
signals  propagating  through  caustic  regions.  The  generalized  DRT 
consists  of  integrating  a  system  of  ordinary  differential 
equations  and  is  based  on  a  natural  extension  of  the  WKBJ  method 
originated  by  Maslov  (Maslov  1972)  in  quantum  mechanics  (see 
Chapman  &  Drummond,  1982,  for  an  alternative  application  of  the 
Maslov-WKBJ  method  in  the  theory  of  seismic  wave  propagation). 

The  wavefield  at  a  receiver  situated  at  some  distance  from  the 
caustics  is  evaluated  by  tracing  the  signals  along  the  rays  con¬ 
necting  the  receiver  to  the  source.  The  field  at  a  point  close  to 
a  caustic  is  given  by  a  diffraction  integral  whose  integrand  can 
be  evaluated  by  tracing  along  a  family  of  rays. 
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AC  a  distance  from  the  caustics  DRT  is  carried  out  by  integrating 
the  usual  system  of  ordinary  differential  equations  along  a  ray. 
The  unknown  functions  include  the  parametric  equation  of  the  ray 
(in  terms  of  the  travel  time  parameter),  slowness,  wavefront  cur¬ 
vature  matrix  and  geometric  ray  spreading  referred  to  the  ray- 
centered  coordinates  as  defined  by  Hanyga  (1982).  In  the  case  of 
S  waves  the  polarization  angle  is  also  included.  This  task  is 
performed  by  a  special  subprogram  of  the  DRT  program. 

Whenever  the  absolute  value  of  the  geometric  ray  spreading  J 
turns  out  to  be  smaller  than  a  positive  number  e  a  canonical 
transformation  in  the  phase  space  of  position  and  slowness  vec¬ 
tors  is  performed  and  a  transformed  version  of  the  DRT  equations 
is  integrated.  This  task  is  performed  by  another  DRT  subprogram. 

As  soon  as  the  geometric  ray  spreading  exceeds  e  >  0  the  original 
DRT  is  resumed. 

The  transformation  of  the  phase  of  the  signal  implicitly  contains 
the  phase  shifts  at  the  caustics.  This  is  a  nontrivial  advantage 
of  the  method  in  the  case  of  anisotropic  media  where  the  phase 
shift  cannot  be  guessed  by  a  simple  argument  (see  Chapman  &  Drum¬ 
mond,  1982,  in  the  case  of  isotropy). 

The  transformations  referred  to  above  are  purely  algebraic.  The 
DRT  system  for  the  caustic  regions  differs  from  the  usual  one 
merely  by  the  number  of  equations  (12  Instead  of  10).  In  the  case 
of  a  ray  which  does  not  exactly  touch  a  caustic  the  transformation 
improves  the  accuracy. 

Boundary  conditions  are  taken  from  Hanyga  (1984a)  and  initial 
conditions  (at  point  sources)  from  Hanyga  (1984d).  DRT  programs 
have  been  designed,  but  testing  has  not  been  completed  yet. 
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Ed£e_dl_f  frjiction_ 

The  DRT  programs  describe..  above  allow  dynamic  ray  tracing  of 
edge-diffracted  rays  from  the  source  to  a  receiver  lying  in  the 
shadow  of  WKBJ  rays  of  a  specified  kind.  An  edge-diffracted  ray 
is  traced  by  the  method  described  in  the  previous  section  from 

the  edge  point  to  the  receiver.  At  the  edge  an  appropriate  ver¬ 
sion  of  Snell's  law  is  satisfied  and  the  amplitude  is  multiplied 
by  the  appropriate  diffraction  coefficient  (Hanyga,  1985). 

For  signal  paths  in  the  transition  zone  of  a  shadow  boundary  uni¬ 
form  asymptotic  formulae  involving  Fresnel  integrals  are  used 
(Achenbach  et  al,  1982;  Hanyga,  1985). 

AnlLsot_ro£y_and_p_rest  res8 

The  references  quoted  in  the  two  previous  sections  deal  with 
arbitrary  anisotropic  and  possibly  prestressed  linear  elastic 
media.  FORTRAN  programs  based  on  these  references  cover  the  cases 
of  isotropic  blocks,  transversely  isotropic  blocks  with  fixed 
orientation  and  with  variable  orientation  of  isotropy  axes. 
Detailed  formulae  for  transversely  isotropic  media  are  taken  from 
Hanyga,  unpublished. 

Klnematix  ra£  tracing  in_3jD 

Preliminary  to  DRT  a  ray  with  a  prescribed  ray  code  is  traced 
from  the  point  source  to  the  receiver.  The  associated  two-point 
boundary  value  problems  are  solved  in  the  following  steps. 

Step  1.  A  large  collection  of  piecewise  straight  rays  is  traced 
through  a  number  of  more  important  interfaces.  These  include  all 
the  efficient  reflectors  in  the  model.  At  each  of  these  inter¬ 
faces  the  ray  splits  into  transmitted,  equal-angle  reflected  and 
two  converted  reflected  rays  (the  latter  are  not  traced  by  some 
programs).  The  rays  are  intercepted  at  their  intersection  with  a 
surface,  a  sphere  or  a  cylinder  and  stored  on  files. 
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Step  2.  The  rays  from  the  above  files  are  used  as  approximate 
rays  in  iterative  programs  tracing  rays  through  homogeneous 

isotropic  blocks.  Polak's  algorithm  (Polak,  1974)  is  used.  Edge- 
diffracted  rays  are  traced  by  the  same  program. 

Step  3.  The  result  of  the  iterative  procedure  of  Step  2  is  fed 
into  a  program  tracing  rays  through  inhomogeneous  and  possibly 
anisotropic  blocks.  The  algorithm  is  based  on  the  idea  of 
parallel  shooting  and  tracing  in  fictitious  time  (Hanyga,  1984a) 
and  uses  Polak's  version  of  the  method  of  secants  (Polak  cit.). 

The  iterative  program  through  isotropic  homogeneous  blocks  takes 
20-30  seconds  in  the  SATM2  models.  The  programs  of  Step  3  take  a 
few  minutes  to  run. 


A.  Hanyga,  NTNF  Postdoc¬ 
torate  Fellow 
(On  leave  of  absence  from 
Institute  of  Geophyscis, 
Warsaw,  Poland) 
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VII. 4 _ Regional  arrays  and  optimum  data  processing  schemes 

In  a  series  of  contributions  (see  reference  list)  we  have 
addressed  various  aspects  of  array  design,  data  processing  sche¬ 
mes  and  past,  present  and  likely  future  developments  in  this 
field.  The  research  to  be  reported  here  is  tied  to  an  evaluation 
of  the  prototype  regional  array  NORESS  (configuration  shown  in 
Fig.  VII.4.1).  In  this  respect  a  two-fold  analyzing  strategy  was 
chosen,  namely,  first  to  estimate  from  NORESS  data  standard  noise 
and  signal  characteristics  like  correlation  functions  and  spectra 
under  a  variety  of  conditions/events  and  secondly,  with  a  basis 
in  Wiener  filtering  theory,  to  evaluate  the  relative  merits  of 
various  array  data  processing  schemes  for  SNR-enchancements .  In 
the  following  the  essence  of  results  obtained  will  be  presented. 

N.°As_e  and  £i£nal_characterj.sdcs 

SNR-spectra  for  events  in  different  tectonic  environments  are 
presented  in  Fig.  VII.4.2a,b,c;  evidently  only  P-waves  with  a 
predominant  shield  path  exhibit  a  significant  amount  of  high- 
frequency  energy  so  as  to  produce  peaks  in  the  SNR-spectra  above 
5  Hz.  Signal  paths  in  oceanic  and  mobile  tectonic  belt  types  of 
lithosphere  appear  to  be  depleted  in  high-frequency  signal 
energy;  the  spectra  here  are  not  too  different  from  those  typical 
of  teleseismic  events. 

Noise  spectra  between  2-10  Hz  appear  to  decay  like  u>“^,  and  be¬ 
sides  that  are  generally  independent  of  time  of  day,  day  of  week 
and  seasonal  variations  as  well.  In  other  words,  the  high- 
frequency  part  of  the  noise  field  is  not  dominated  by  cultural 
and  microseisraic  noise  contributions. 

Noise  correlations  as  functions  of  sensor  separations  have  to  be 
calculated  under  a  variety  of  conditions,  that  is,  sensor  time 
lags  tied  to  both  vertical  and  horizontal  beams.  Despite  large 
fluctuations  in  individual  estimates  here,  significant  differen¬ 
ces  exist  as  a  function  of  beam  location,  e.g.,  see  Figs. 

VII. 4. 3a  and  b.  This  feature  is  rather  obvious  from  Table 
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VII. 4.1,  where  SNR-enhancements  as  a  function  of  phase  velocity 
and  azimuth  for  different  bandpass  filters  and  processing  schemes 
amounts  to  approximately  3  dB! 

Arra^  da_ta_process ing_  schemes. 

Within  the  general  class  of  Wiener  filtering  3  array  data  pro¬ 
cessing  schemes  have  been  considered,  namely,  standard  beam¬ 
forming,  optimum  weighting  and  maximum  likelihood  (ML)  filtering. 
The  mathematics  of  these  approaches  are  detailed  in  Ingate  et  al 
(1985).  Standard  beamforming  is  optimum  for  uncorrelated  noise 
(/N  gain)  which  in  general  takes  place  above  4  Hz  for  the  NORESS 
array  configuration.  Optimum  weighting  may  give  an  additional 
2-4  dB  gain  for  correlated  noise,  while  ML-processing  is  clearly 
superior  vis-d-vis  the  other  two  processing  schemes,  as  evident 
from  Fig.  VII. 4. 4.  The  well-known  problem  with  more  advanced 
array  data  processing  schemes  is  their  rather  severe  computer 
loads.  For  example,  for  the  optimum  weighting  scheme  the  noise 
covariance  matrix  has  to  be  updated  each  2. 5-5.0  sec  because  of  a 
general  noise  non-stationarity  as  illustrated  in  Fig.  VII. 4. 5. 

Discussion 

After  extensive  analysis  and  processing  of  the  NORESS  data,  some 
remarks  on  regional  array  design  and  operation  may  be  justified. 
First  of  all,  array  design  is  problematic  as  the  operational 
bandwidth  reflecting  SNR-spectra  from  different  tectonic  regimes 
should  cover  roughly  the  3-8  Hz  range  in  contrast  to  "tele- 
seismic"  arrays,  where  the  bandwidth  is  1-2  Hz.  To  ensure 
adequate  seismicity  surveillance  capabilities  in  the  3-8  Hz  band 
very  many  beams  have  to  be  formed,  otherwise  signal  decorrelation 
losses  would  be  severe,  even  when  the  data  processing  is 
restricted  to  simple  beamforming.  Indeed,  if  relatively  better 
performance  is  desired,  this  may  most  easily  be  achieved  on  the 
hardware  side  by  adding  more  instruments  and/or  using  the  "analog 
sensor  clustering"  commonly  used  in  seismic  land  surveys.  As 
regards  array  conf iguraiton  per  se,  a  reasonably  flat  sensor 
spacing  distribution  has  to  be  ensured  if  we  want  a  reasonably 
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flat  array  response  as  a  function  of  frequency.  For  example,  from 
Fig.  VII. 4. 4  we  see  that  the  NORESS  performance  is  far  from  opti¬ 
mum  below  4  Hz  due  to  high  noise  correlation  and  likewise  beyond 
6-8  Hz  due  to  poor  signal  correlations. 


Arrays  are  relatively  costly  to  build  and  operate  so  there  have 
been  limited  possibilities  for  practical  experimentation. 

However,  taking  advantage  of  recent  advances  in  microprocessor 
and  communication  technology,  it  is  feasible  to  construct  and 
deploy  small,  inexpensive  arrays  as  discussed  in  detail  by  Huse- 
bye  et  al.  This  would  naturally  add  to  the  flexibility  of  prac¬ 
tical  array  experiments,  and  also  increase  the  number  of  arrays 
being  operative.  Apparently,  most  of  the  seismological  community 
takes  part  in  these  new  developments  as  the  newly  formed  IRIS- 
consortium  plans  for  a  digital,  global  seismograph  network  and 
mobile  seismic  arrays  comprising  hundreds  (thousands)  of 
instruments.  These  and  likely  future  trends  in  digital  seismome- 
try  are  discussed  in  all  papers  referenced  below. 

Result_summa ry  o_f  the;  j>rot ot^pe  TOREj3S_evaluat^_o_n 

The  major  problem  here  was  that  of  suppressing  ambient  noise,  and 

major  results  obtained  are  as  follows. 

Signal  and  noise  field  characteristics 

SNR  spectra  peak  at  ca  2,  3  and  6  Hz,  respectively,  for 
events  whose  signal  paths  are  predominantly  oceanic, 
mobile  cratonic  belt  and  shield.  Teleseismic  arrays  like 
LASA  and  NORSR  operate  generally  in  the  narrow  1-2  Hz 
band. 

Above  2  Hz  noise  spectra  are  essentially  time  invariant 
and  besides  have  approximately  a  decay  rate  in  the 
2-10  Hz  band. 


Noise  correlation  distance  is  clearly  frequency  depen¬ 
dent  and  is  such  that  noise  suppression  by  simple  beam¬ 
forming  will  be  degraded  relative  to  /N  up  to  about 
4  Hz. 

The  noise  field  cannot  be  considered  stationary  even  for 
short  windows  of  2. 5-5.0  sec.  This  is  most  easily  seen 
in  standard  beamforming  gain  variation  up  to  3  dB  as  a 
function  of  beam  location  or  azimuth  and  phase  velocity. 
Signal  correlations  for  P  and  Lg  phases  decrease  with 
increasing  sensor  separation  becoming  ca  0.5  at  about 
3  km.  This  in  combination  with  steering  delay  errors, 
including  insufficient  number  of  beams  deployed,  would 
severely  affect  array  performance  above  6-8  Hz. 
Horizontal  seismometer  recordings  exhibit  noise  and 
signal  correlations,  etc.,  roughly  similar  to  that 
observed  for  vertical  instruments. 

Noise  suppression  schemes  -  variants  of  Wiener  filter  theory 

Simple  beamforraing  is  optimal  or  /W  for  uncorrelated 
noise,  but  this  cannot  be  achieved  using  all  instruments 
due  to  the  relatively  small  NORESS  array  aperture  of 
3  km. 

Optimum  weighting  based  on  characteristics  as  manifested 
in  the  noise  covariance  matrix  (multichannel)  will  give 
an  additional  gain  of  2-4  dB  relative  to  standard  beam¬ 
forraing  up  to  about  5  Hz.  A  simplification  here, 
reducing  these  weights  to  the  0/ 1-type,  gives  1-2  dB 
less  gains  but  still  better  than  standard  beamforming. 
Maximum  likelihood  filtering,  computationally  demanding, 
gave  approximately  /TT  gains  even  for  lower  frequencies 
where  simple  beamforming  was  inefficient. 

On  the  basis  of  the  above  results  and  technical  considerations  as 

well,  our  conclusion  is  that  seismic  surveillance  is  most  effec¬ 
tively  attained  by  deploying  rather  many  seismometers  (regional 

arrays)  combined  with  simple  processing  schemes  rather  than  few 
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Instruments  combined  with  sophisticated  processing  schemes  like 
ML-filtering  -  simply  because  investments  in  array  hardware  are 
considered  more  cost-effective  than  array  software  investments. 


G.S.  Husebye 
S.F.  Ingate  (NTNF  Post¬ 
doctorate  Fellow) 

A.  Christof fersson 
(Uppsala  University) 
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Table  VII. 4.1  Standard  beamforming  and  optimum  weighting  gains 
as  functions  of  primarily  azimuth  but  also  phase 
velocity  and  filter  passband.  The  thirteen 
NORESS-sensors  used  were  nos.  1,  5  to  9  and  10  to 
16.  Azimuthal  gain  variation  amount  to  approxima¬ 
tely  2  dB,  which  may  increase  1  dB  if  finer  azi¬ 
muth  sampling  intervals  of  30  deg  had  been  used. 


Fig.  VII. 4.1  The  experimental  NORESS  array  located  near  NORSAR 


site  06C.  Instruments  at  sites  indicated  by  open 
circles  will  be  in  operation  in  late  1984. 


NORMALIZED  POWER  (dB)  NORMALIZED  POWER  (dB) 


Fig.  VII. 4. 3  Mean  normalized  correlation  curves  for  noise  on  the 


NORESS  array  as  a  function  of  inter-sensor  spacing  and 
frequency,  a)  vertical  beam  lags;  b)  horizontal  beam 
lags  (azimuth  -  0°,  velocity  -  4.2  km/s).  Each  curve  is 
an  average  of  100  sec  windows  sampled  every  5  hours  for 
two  days. 
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Fig.  VII. 4. 5  Display  of  standard  and  optimum  beam  traces  formed  during  the 
processing  of  noise  recorded  by  NORESS.  The  corresponding 
gain  function  is  not  smooth,  but  rather  fluctuates  rapidly  in 
certain  time  intervals  which  reflects  a  degree  of  non- 
stationarity  in  the  ambient  noise  field.  The  first  10  seconds 
of  data  was  used  to  stabilize  the  covariance  matrix,  and 
consequently,  not  included  in  the  processing.  Amplitude 
scales  for  the  beam  traces  are  in  amplitude  units,  while  the 
amplitude  scale  for  the  gain  function  is  in  dB. 
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VII. 5  High-frequency  P  wave  attenuation  from  Central  Asia  to 
Norway 

It  has  long  been  known  that  the  P  waves  recorded  in  northern  Europe 
from  earthquakes  and  explosions  in  Eurasia  are  often  quite  rich  in 
high-frequency  energy.  A  detailed  quantification  of  signal  trans¬ 
mission  characteristics  is  then  possible  by  studying  the  spectral 
characteristics  over  a  broad  frequency  band.  In  this  study  the  spectra 
of  NORSAR  recordings  of  some  E.  Kazakh  events  are  used  to  develop  a 
model  for  the  average  Q  along  this  path.  Since  the  epicentral  distance 
is  about  4200  km,  this  average  Q  is  for  paths  that  penetrate  the  earth 
to  a  depth  of  about  900  km.  The  data  were  processed  with  the  technique 
applied  by  Bache  et  al  (1984).  For  each  seismogram  energy  density 
spectra  are  computed  for  very  short  (typically  2.2  to  2.5  seconds) 
time  windows  Isolating  the  first  arriving  compressional  wave,  and  the 
power  spectrum  of  a  noise  window  just  before  the  signal  is  subtracted. 
Event  spectra  are  then  obtained  by  averaging  over  several  elements  of 
the  array  as  well  as  similar  events  in  a  small  source  region,  and  the 
spectra  obtained  this  way  are  quite  smooth  and  are  almost  entirely 
shaped  by  the  average  source  spectrum  and  the  effect  of  attenuation 
along  the  path.  Accurate  attenuation  estimates  depend  on  accurate 
corrections  for  the  source,  assuming  that  the  source  spectrum  for  the 
direct  P  wave  should  be  roughly  constant  up  to  some  corner  frequency, 
above  which  it  is  proportional  to  f~n.  For  explosions,  most  of  the 
evidence  supports  n  ■  2,  and  the  NORSAR  results  appear  to  confirm  that 
value. 

The  spectra  computed  from  five  8  Hz  seismometers  for  two  large  Degelen 
events  are  shown  in  Fig.  VII. 5.1.  Also  shown  is  the  noise  spectrum 
computed  from  the  average  noise  power  (based  on  a  3.8  second  sample  of 
noise  before  the  signal  on  each  element).  The  signal  spectrum  is 
essentially  flat  from  3  to  beyond  7  Hz  when  plotted  with  the  f^  source 
correction,  and  the  spectra  for  the  other  events  look  much  the  same. 


This  is  an  extraordinary  result,  as  the  attenuation  has  almost  no 
frequency-dependent  effect  between  3  and  7  Hz.  Since  attenuation  can¬ 
not  be  negative,  spectra  like  this  confirm  the  assumption  that  the 
source  is  proportional  to  f-^  above  the  corner  frequency. 

In  Fig.  VII. V. 2  we  show  the  spectra  for  6  other  events.  Four  are 
smaller  than  the  Degelen  explosions  and  two  are  explosions  in  other 
paths  of  the  E.  Kazakh  test  site.  The  spectra  for  the  Degelen  events 
are  all  similar  to  those  in  Fig.  VII. V. 2,  except  that  they  have  the 
different  low-frequency  behavior  expected  for  smaller  (i.e.,  higher 
corner  frequency)  events.  The  maximum  signal/noise  occurs  around  7  Hz 
and  decreases  at  higher  frequencies.  The  events  in  the  other  two  areas 
seem  to  depart  from  the  pattern  in  that  the  signal  spectrum  does  not 
decay  so  obviously  above  7-8  Hz. 

Three  major  features  of  these  spectra  determine  the  attenuation  model 
for  the  E.  Kazakh-NORSAR  path.  First,  at  low  frequencies  they  are 
influenced  by  source  corner  frequency  effects  and  a  correction  for 
them  must  be  introduced.  Second,  the  source-corrected  spectra  are 
essentially  flat  in  the  band  between  3  and  7  Hz.  Finally,  the  spectra 
decrease  above  7-8  Hz,  and  we  see  in  Fig.  VII. 5. 3  that  the  line  repre¬ 
senting  t*  =  0.14  is  a  reasonable  approximation  for  the  rate  of  decay. 

For  frequencies  up  to  7  Hz,  the  attenuation  effects  can  be  represented 
by  an  absorption  band  model.  This  is  best  done  with  a  "path  average" 
spectrum  obtained  by  averaging  several  spectra  (here  six),  as  shown  in 
Fig.  VII. 3. 4,  obtaining  a  spectrum  which  essentially  is  a  smoothed 
version  of  the  single  event  spectra  in  Fig.  VII.5.1.  In  Fig.  VII. 5. 4 
we  also  show  an  estimate  for  the  spectrum  after  correcting  for  source 
corner  frequency  effects  and  the  spectra  for  several  attenuation 
models.  The  corner  frequency  correction  is  based  on  spectra  computed 
by  Bache  et  al  (1984),  utilizing  the  spectral  difference  between  Dege- 
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len  and  low  corner  frequency  Shagan  events  to  approximately  correct 
the  Degelen  spectra.  A  roughly  "corrected"  spectrum  (Fig.  VII. 5. 4) 
shows  that  the  slope  changes  below  3  Hz,  even  though  this  "corrected" 
spectrum  still  underestimates  attenuation  effects  at  low  frequency. 

In  Fig.  VII. 5. 4  the  "corrected"  NORSAR  spectrum  is  plotted  with  the 
spectral  effect  of  several  absorption  band  models  (Minster,  1978).  For 
this  application  the  model  has  two  parameters,  tQ*,  which  is  the  ratio 
of  travel  time  to  average  path  Q  at  long  periods,  and  Tm,  which  speci¬ 
fies  where  the  Q-^  begins  to  decrease  with  frequency.  The  key  spectral 
features  are  the  frequency  dependence  below  3  Hz  and  the  nearly  hori¬ 
zontal  slope  up  to  7  Hz.  The  best  fit  seems  to  be  with  models  having 
t0*  =  0.6  seconds  (not  much  constrained  by  these  data)  and  Tm  slightly 
greater  than  0.05  seconds.  For  such  models  the  effect  of  attenuation 
above  3  Hz  reduces  to  multiplication  of  the  spectrum  by  a  constant, 
which  is  about  0.04  for  tQ*  *  0.6,  tm  *  0.06;  and  0.07  for  tD*  =  0.5, 
tm  -  0.06.  Obviously,  the  value  of  this  constant  depends  on  the  choice 
of  t0*,  characterized  by  the  attenuation  at  long  periods. 

An  absorption  band  model  with  0.05  <  <  0.1  and  t0*  =  0.6  was  found 

by  Bache  et  ai  (1984)  to  represent  the  low  frequency  attenuation  on 
the  paths  from  eastern  Kazakhstan  to  four  UKAEA  arrays,  with  frequen¬ 
cies  above  3  Hz,  requiring  a  second  Q-^  to  explain  the  constant  rate 
of  spectral  decay  seen  at  those  frequencies.  In  view  of  Figs.  VII. 5. 2 
VII. 5. 3,  the  NORSAR  path  data  also  require  a  second  Q  superimposed  on 
the  absorption  band  model,  but  its  effect  is  not  felt  until  7  Hz  and 
beyond.  Thus,  similar  two-part  Q  models  seem  to  fit  the  data  for  all 
five  paths  (with  minor  adjustments  in  the  t0*  and  xm),  with  the  major 
difference  being  that  the  effect  of  the  second  Q  is  not  felt  until 
higher  frequencies  at  NORSAR. 

A  key  feature  of  the  NORSAR  spectra  is  thus  that  the  effect  of  attenu¬ 
ation  is  little  more  than  multiplication  by  a  constant  for  frequen- 


cies  between  3  and  7  Hz.  Physical  mechanisms  for  absorption  must  be 
characterized  by  Q  =  f  above  some  high-frequency  limit  on  the  relaxa¬ 
tion  spectrum,  so  absorption  band  models  (e.g.,  Anderson  and  Given, 
1982)  anticipate  that  attenuation  must  reduce  to  a  constant  multiplier 
above  some  limiting  frequency.  These  NORSAR  spectra  provide  an  unam¬ 
biguous  observation  that  such  behavior  does  indeed  exist.  For  most 
paths  it  is  difficult  to  observe  apparently  because  attenuation  by 
scattering  masks  it. 

It  is  not  known  whether  or  not  the  E.  Kazakh-NORSAR  path  is  typical  of 
the  Eurasian  continent.  If  it  is,  the  detection  capability  of  systems 
to  monitor  nuclear  explosion  testing  in  this  region  could  be  much 
improved  by  greater  emphasis  on  high-frequency  data.  An  estimate  of 
the  anticipated  improvement  was  recently  made  by  Evernden  et  al 
(1984),  but  it  should  be  noted  that  these  authors  made  much  more  opti¬ 
mistic  assumptions  about  the  signal/noise  at  high  frequencies  than  can 
be  supported  by  the  NORSAR  data  presented  here.  Some  improvement  of 
the  NORSAR  signal/noise  at  high  frequencies  may  be  possible  by  using 
arrays  designed  specifically  for  this  purpose  (e.g.,  Mykkeltveit  et 
al,  1983),  but  the  scattering  attenuation  may  be  impossible  to  over¬ 
come  since  it  is  probably  due  to  inhomogeneous  structure  throughout 
the  lithosphere.  Regional  phases  like  Pn  also  seem  likely  to  be  at¬ 
tenuated  by  scattering  at  high  frequencies.  While  better  use  of  high- 
frequency  data  may  improve  detection,  this  is  far  easier  than 
identifying  an  event,  and  more  work  remains  to  be  done  until  one  knows 
how  much  the  capability  for  improving  event  identification  can  be 
Improved  by  exploiting  these  high-frequency  data. 
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More  details  about  the  research  reported  on  here  are  given  in  Bache 
and  Bungum  (1984). 

T.C.  Bache  (Science  Appli¬ 
cations  International, 

La  Jolla) 

H .  Bungum 
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Fig.  VII. 5.1  The  spectra  for  two  of  the  events  of  (Table  2  ?)  are 

shown  together  with  the  average  noise  window  just  before 
the  P  signal.  A  horizontal  dashed  line  is  sketched 
through  each  signal  spectrum  and  the  NORSAR  noise  model 
of  Bungum  (1983)  is  sketched  through  the  noise. 
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Fig.  VII. 5. 3  Signal  and  noise  spectra  are  plotted  for  an  event  digi¬ 
tized  at  two  different  sampling  rates.  The  higher 
sampling  rate  spectra  are  then  replotted  on  a  different 
frequency  scale  (bottom).  On  the  latter  a  line  with  a 
slope  corresponding  to  t*  ■  0.14  is  drawn  through  the 
high-frequency  portion  of  the  signal  spectrum. 


2.0 


4.0  6.0 

FREQUENCY,  Hz 


80 


10.0 


Fig.  VII. 5. 4  The  bottom  sketch  shows  the  path  average  spectrum 
together  with  a  horizontal  line  and  a  dotted  line 
showing  this  spectrum  corrected  for  corner  frequency 
effects.  The  top  two  sketches  show  this  corrected 
spectrum  with  the  spectral  effect  of  several  absorption 
band  models  for  attenuation. 


VII. 6  Teleselsmic  detection  at  high  frequencies  using  NORSAR  data 

This  paper  summarizes  initial  results  from  NORSAR  noise  and  signal 
studies  aimed  at  identifying  possibilities  to  improve  the  detection 
capability  of  the  NORSAR  array.  Particular  emphasis  is  given  to  high- 
frequency  signals  (2-5  Hz),  which  are  typical  for  NORSAR  recordings  of 
many  Eurasian  events. 

Noise  characteristics .  Bungura  and  Mykkeltveit  (1984)  have  summarized 
the  main  features  of  seismic  noise  spectra  observed  at  NORSAR.  These 
spectra  are  characterized  by  a  very  strong  microseismic  peak  at 
0.2-0. 3  Hz,  especially  during  North  Sea  storm  activity.  The  spectral 
slope  is  very  steep  (50  dB/decade)  up  to  at  least  40  Hz,  and  the 
NORSAR  noise  at  high  frequencies  is  among  the  lowest  observed.  Above  2 
Hz,  there  is  little  difference  between  "high  noise"  and  "low  noise" 
conditions,  except  that  cultural  activities  cause  increases  in  high 
frequency  noise  during  day  time.  Because  of  the  limited  resolution  and 
gain  ranging  applied  at  NORSAR,  the  spectra  actually  observed  at 
NORSAR  instruments  are  biased  high  above  3  Hz  as  discussed  by  Bungum 
(1983). 

Signal  characteristics.  Main  features,  discussed  by  Ringdal  and 
Husebye  (1982)  are  a)  large  and  region-dependent  variation  in  signal 
levels  across  NORSAR,  often  spanning  an  order  of  magnitude  across  the 
array,  b)  significant  deviations  from  a  plane  wavefront  model,  making 
steering  delay  corrections  necessary  in  the  beamforming  process,  and 
c)  significant  energy  at  high  frequencies  (>  2  Hz)  from  Eurasia,  espe¬ 
cially  for  underground  explosions. 

Detection  algorithms .  Current  Detection  Processor  (DP)  procedure  at 
NORSAR  (Ringdal,  1981)  is  to  form  conventional  array  beams  with  a 
filter  of  1.2-3. 2  Hz  and  envelope  beams  with  a  filter  of  1.6-3. 2  Hz. 
These  beams  are  then  subjected  to  a  linear  STA/LTA  detector  and  passed 
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through  a  grouping  algorithm  to  delete  side  lobe  detections.  In  the 
present  study,  we  have  analyzed  SNR  of  a  number  of  presumed  explosions 
from  Eurasia  by  applying  8  predetermined  bandpass  filters  to: 

(a)  An  array  beam  using  current  DP  steering  delays  (0.1  sec 
accuracy) 

(b)  An  array  beam  using  recomputed  delays  (0.05  sec  accuracy) 

(c)  An  optimum  weighted  array  beam,  with  steering  delays  as  in  (b) 
and  each  channel  weighted  by  S/N^  (S  *  signal  amplitude,  N  * 
noise  amplitude) 

(d)  The  best  subarray  beam  for  the  particular  event. 

Examples  of  results  are  shown  for  two  events  in  Figs.  VII. 6.1  and 
VII. 6. 2.  These  results,  which  are  largely  confirmed  (with  few  excep¬ 
tions)  by  the  other  events  analyzed,  give  that: 

(a)  Significant  gains  are  possible  by  applying  more  high-frequent 
filters  than  currently  done. 

(b)  The  weighted  array  beam  is  generally  best,  while  the  best 
subarray  beam  SNR  is  quite  close  (within  2-3  dB  average). 

(c)  Gains  compared  to  current  processing  can  reach  10-20  dB  for 
high-frequent  signals,  and  somewhat  less  for  signals  of  low 
dominant  frequency. 

Detection  performance  for  presumed  explosions.  By  comparing  the  SNR  of 
the  best  beam  to  the  event  magnitude,  it  is  possible  to  estimate 
approximately  the  optimum  N0RSAR  thresholds,  i.e.,  thresholds  that 
could  be  achieved  if  processing  as  indicated  above  were  implemented. 

We  have  found  that  (cf.  Figs.  VII. 6. 3  and  VII. 6. 4)  these 
"instantaneous"  detection  thresholds  are: 
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West  of  Ural  Mountains:  mb  2. 0-2. 5  (possibly  better) 

Caspian  Area  2. 0-2. 5 

Semipalatinsk  2. 5-3.0 

Siberia  mb  2. 5-3. 5 

It  is  noteworthy  that  there  is  a  large  range  in  thresholds  across 
Eurasia,  and  this  can  be  interpreted  as  being  due  to  "source  focusing" 
effects  similar  to  the  "receiver  focusing"  effects  routinely  observed 
at  NORSAR. 

Detection  in  the  coda  of  an  earthquake.  An  example  of  the  importance 
of  high-frequency  detection  in  this  connection  is  given  in  Fig.  VII. 6. 5, 
where  NORSAR  06C06  signals  from  an  mb  5.8  earthquake  near  Kamchatka 
are  followed  one  minute  later  by  signals  from  a  presumed  explosion 
(mb  3.8)  at  Semipalatinsk.  With  the  standard  1.2-3. 2  Hz  filter,  the 
latter  signal  is  completely  masked  by  the  coda,  whereas  a  filter  of 
3. 2-5. 2  Hz  shows  the  explosion  signal  dominating  that  of  the 
earthquake.  This  figure  is  based  on  a  single  instrument,  and  further 
improvements  are  of  course  possible  by  multichannel  processing. 

In  conclusion,  the  demonstrated  gains  from  high-frequency  processing 
of  teleseismic  explosions  clearly  warrant  a  modification  in  the  pre¬ 
sent  NORSAR  detection  algorithms.  The  practical  implementation  of 
these  changes  is  now  the  subject  of  further  study. 

F.  Ringdal 
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,  VII. 6.1  Observed  SNR  at  NORSAR  for  different  filters  and  bea 
forming  methods,  as  explained  In  the  text.  The  case 
Is  an  mh  A. 5  Semipalatinsk  presumed  explosion. 
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Fig. VII. 6. 2  Same  as  Fig.  VII. 6.1,  but  for  a  presumed  explosion  from 
Che  Azgir  area. 
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Fig.  VII. 6. 3  NORSAR  array  beam  (top)  and  a  scaled  beam  added  to  pre¬ 
ceding  noise,  and  indicating  the  "optimum"  instantaneous 
dectection  threshold.  This  figure  corresponds  to  the  event 
of  Fig.  VII. 6.1.  Note  that  this  is  a  "best  case"  event, 
and  that  typical  scaled  thresholds  for  other  Semipala- 
tinsk  events  range  from  m^  2.5  to  3.0. 
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Fig.  VII. 6. 5 
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Example  of  NORSAR  recordings  instrument  06C06  from  a  small 
presumed  explosion  with  signal  arrival  one  minute  later 
than  that  of  a  preceding  large  earthquake.  Note  that 
the  explosion  signal  is  not  visible  on  the  top  trace 
(1.2-3. 2  Hz  filter),  but  can  be  clearly  seen  on  the 
bottom  trace  (3. 2-5. 2  Hz  filter). 
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VII. 7.1  NORESS  Regional  Array  developments 

Evaluation  of  the  performance  of  the  newly  established  NORESS  array  in 
Norway  has  been  initiated.  In  particular,  real  time  processing 
algorithms  aimed  at  detection  and  location  of  events  at  regional 
distances  have  been  developed  and  tested. 

The  design  criteria  and  associated  specifications  on  the  performance 
of  the  new  array  were  originally  formulated  as  follows: 

Detection:  The  array  should  provide  optimum  SNR  gain  by 
beamforming. 

Location:  The  array  should  have  a  narrow  main  lobe  and  small  side 
lobes. 

Performance  as  a  function  of  signal  frequency:  The  array  perfor¬ 
mance  should  be  good  or  at  least  reasonable  for  a  wide  range  of 
frequencies,  typical  of  regional  wave  propagation. 

-  Research  aspects:  The  array  should  have  an  optimum  coarray 

pattern.  Add  as  many  horizontal  component  seismometers  as  can  be 
afforded  (within  the  geometry  of  the  vertical  array). 

The  array  geometry  of  Fig. VII. 7.1  was  worked  out  in  close  cooperation 
between  Sandla,  LLNL  and  NORSAR.  It  represents  the  best  compromise  we 
found  between  the  partly  conflicting  demands  made  above  on  array 
performance.  In  the  design  work,  we  were  able  to  draw  heavily  on  the 
experience  gained  from  several  provisional  array  installations  in 
Norway.  For  example,  actually  observed  correlations  for  signals  and 
noise  were  used  to  optimize  the  SNR  gain.  For  evaluation  of  the  loca¬ 
tion  capability  of  the  new  array,  beam  patterns  were  computed  for 
realistic  signal  correlations.  The  geometry  of  Fig.  VII. 7.1  contains 
useful  subgeometries,  which  are  close  to  optimum  for  different  fre¬ 
quencies  within  the  range  2-10  Hz. 

In  order  to  have  an  experimental  check  on  the  proposed  geometry,  a 
provisional  21-element  array  was  installed  during  the  summer  of  1983. 
Analysis  of  data  from  this  installation  has  confirmed  the  correlation 
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curves  for  signals  and  noise  on  which  the  design  work  was  based. 

Testing  of  the  RONAPP  processing  package  (Mykkeltveit  &  Bungura,  1984) 
for  regional  events  recorded  on  this  array,  shows  that  the  arrival 
azimuths  can  be  estimated  to  within  a  few  degrees  of  the  true  value. 

Results  from  the  1983  provisional  installation  basically  confirmed  the 
projected  capabilities  of  the  new  array,  and  it  was  finally  decided  to 
deploy  an  array  with  the  geometry  given  in  Fig.  VII.7.1. 

The  final  array  installed  during  the  summer  of  1984  comprises  25  short 
period  elements,  4  out  of  which  are  3-axis  deployments  with  the 
remaining  21  elements  consisting  of  vertical-motion  seismometers  only. 

In  addition,  there  is  a  broad-band  3-axis  system.  Sampling  rates  are 
40  Hz  for  the  short  period  channels  and  10  Hz  (intermediate  period) 
and  1  Hz  (long  period)  for  the  broadband  system.  Data  from  the  new 
array  are  transmitted  to  Kjeller  via  a  leased  land  line  and  via 
satellite  to  the  U.S.  Preliminary  analysis  of  the  new  data  has  given 
promising  results.  Fig.  VII. 7. 2  gives  an  example  of  the  automatic  output 
from  the  on-line  event  processing  of  data  from  the  new  array. 

Several  improvements  to  the  on-line  processing  package  are  under  way. 

For  instance,  the  success  of  the  location  part  of  the  RONAPP  pro¬ 
cessing  package  depends  on  a  proper  identification  of  the  secondary 
phase  associated  with  the  P  arrival.  This  secondary  phase  is  either  Lg 
or  Sn  and  a  priori  knowledge  on  the  occurrence  of  these  phases  on  the 
NORESS  records  as  a  function  of  source  region  will  be  incorporated  in 
RONAPP.  For  this  end,  regional  data  recorded  at  NORSAR  over  the  past 
13  years  have  been  examined.  The  general  picture  is  that  Lg  is  the 
dominant  secondary  arrival,  but  notable  exceptions  have  been  iden¬ 
tified,  corresponding  to  source  regions  for  which  Sn  is  the  more  pro¬ 
minent  or  even  sole  secondary  phase.  Fig.  VII. 7. 3  gives  an  example  of 
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a  tectonic  feature  in  the  North  Sea  that  is  an  effective  barrier  to  Lg 
wave  propagation  (Kennett  &  Mykkeltveit,  198A). 


S.  Mykkeltveit 
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Fig.  VII. 7. 2a  Vertical  component  records  from  the  new  NORESS  array 
for  a  local  event  detected  and  located  by  RONAPP.  The 
two  arrivals  detected  and  associated  are  marked  by 
arrows.  Bottom:  Unfiltered  and  filtered  beams  for  the 
two  detections. 
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Fig.  VII. 7. 2b  Results  of  frequency-wavenumber  analysis  of  the  indi¬ 
cated  time  windows  around  the  two  detections  in  Fig. 
VII. 7. 2a. 


Fig.  VII. 7. 3  Map  of  the  central  North  Sea  Basin  (top)  and  main  struc¬ 
tural  features  (middle).  NORSAR  records  from  shots 
located  at  positions  indicated  in  the  map.  Arrows  indi¬ 
cate  a  group  velocity  of  3.5  kins'"  *  corresponding  to  Lg 
wave  arrivals. 


